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The cosmic microwave background provides constraints on the annihilation and decay of light 
dark matter at redshifts between 100 and 1000, the strength of which depends upon the fraction 
of energy ending up in the form of electrons and photons. The resulting constraints are usually 
presented for a limited selection of annihilation and decay channels. Here we provide constraints on 
the annihilation cross section and decay rate, at discrete values of the dark matter mass m x , for all 
the annihilation and decay channels whose secondary spectra have been computed using PYTHIA 
in arXiv:1012.4515 ("PPPC4DMID: A Poor Particle Physicist Cookbook for Dark Matter Indirect 
Detection"), namely e, fx, t, V — ¥ e, V — > fi, V — ► t, u, d s, c, b, t, 7, g, W, Z and h. By 
interpolating in mass, these can be used to find the CMB constraints and likelihood functions from 
WMAP7 and Planck for a wide range of dark matter models, including those with annihilation or 
decay into a linear combination of different channels. 



The temperature and polarization fluctuations of the 
cosmic microwave background (CMB) are well known to 
be sensitive to the redshift of recombination, z w 1100, 
as this determines the surface of last scattering. If dark 
matter annihilation or decay deposits electromagnetic en- 
ergy in the primordial plasma after z ~ 1100, it can delay 
recombination and/or contribute to reionization, lead- 
ing to distortions in the CMB [1-19]. This is especially 
constraining for light dark matter x with mass m x < 
lOGeV (as its number density is greater than that of 
heavier dark matter) , and if annihilation or decay is into 
electrons; in that case the current limit on the annihi- 
lation cross section is close to the standard relic density 
value (av) = 3 x 10~ 26 cm 3 /s. Apart from the special 
case of monochromatic photons, constraints are generi- 
cally weaker for annihilations or decays into other par- 
ticles, as those channels all involve substantial eventual 
yields into neutrinos and/or hadrons, neither of which 
efficiently transfer energy to the primordial gas (neutri- 
nos because they are weakly-interacting, hadrons because 
they are strongly penetrating [4]). 

Recent progress has been made in refining and sys- 
tematizing the CMB bounds on dark matter annihilation 
and decay [14, 17]. A key quantity for determining the 
constraint on a given model is the efficiency f(z) for pro- 
ducing ionizing radiation, as a function of redshift z. For 
annihilations, f(z) is denned in terms of the electromag- 
netic power injected per unit volume, 



dE 
dt dV 



= /(z)^l£ M p2c»(l + *) a 



(1) 



where p c is the critical mass density of the universe to- 
day and f^DM is the fraction in dark matter. For perfect 
efficiency, / = 1, this is just twice the the DM mass times 
the annihilation rate per unit volume. Note that no fac- 
tors of two remain in (1), as the annihilation rate itself 



contains a factor of one half for annihilation of identical 
particles [20], which cancels the factor of two from the 
release of twice the DM rest mass in each annihilation 
event. Here we assume Majorana DM; were the DM not 
its own antiparticle, (1) would need to be divided by a 
further factor of two, and IJdm interpreted as the total 
mass fraction of DM + anti-DM particles. 

Ref. [17] has provided transfer functions Ti(z',z,E) 
that determine the contribution to f(z) from particles i = 
7 or pairs injected at redshift z' with energy E. The 
injected 7 or may be primary products of dark matter 
annihilation or decay, or they may be secondary particles 
arising from the showering and decay of the primary ones. 
If fa — dNi/dE is the electron/photon spectrum of state 
i, normalized such that / dE E fa is the fraction of initial 
dark matter mass converted into energy in that state, 
then 2 



f(z)= j dz' J dEEfi(z',z,E)fa{E). 



(2) 



The departure from perfect efficiency, / = 1, is due to 
the fraction of initial energy that ends up in neutrinos 
(or possibly more exotic invisible final states, which we 
do not explicitly consider here). 

To find the spectra fa, Monte Carlo computations us- 
ing an event generator such as PYTHIA [21, 22] or HER- 
WIG [23] are required. These simulations are numeri- 
cally intensive, but the results have been carried out and 
made available for a range of primary annihilation and 
decay channels in ref. [24]. The spectra obtained from 
PYTHIA 8.135 and HERWIG 6.510 generally agree to 
within about 20% for most primary annihilation/decay 
channels, except for g (gluons), where the choice of event 
generator results in up to a factor of two uncertainty in 
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TABLE I: / 8 g values as a function of WIMP mass m x (in 
GeV) and primary annihilation channel (for example "e" de- 
notes XX ~^ e 6, while "V e" denotes \\ — > VV , followed 
by V — > ee), for computing WMAP7 (left) and projected 
Planck (right) constraints. Null entries indicate that m x is 
below threshold. Bottom row covers extra range of masses 
where \X ~~ >" 11 constraints undergo a transition in fig. 2. 
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FIG. 1: First three principal components for annihilating dark 
matter, from http://nebel.rc.fas.harvard.edu/epsilon/. Solid 
curves are for WMAP7 and dashed for anticipated Planck 
data. 



electron and photon yields. At lower energies PYTHIA 
also gives substantially larger photon yields for leptonic 
and gauge boson channels, a consequence of the omis- 
sion of (QED) lepton final state radiation and fermion 
pair production in HERWIG. We therefore exclusively 
employ the PYTHIA results of rcf. [24], which have the 
additional benefit of also including electroweak correc- 
tions [25]. 

The purpose of the present note is to provide a set of 



results that can be easily interpolated to find CMB con- 
straints and likelihood functions for dark matter models 
that annihilate or decay to arbitrary final states with a 
constant cross-section or decay rate. We first carry out 
the computation of (2) for a full range of dark matter 
masses and annihilation and decay channels, and then 
translate each result into a single number that encodes 
the energy injection history for that mass and final state. 
For annihilation we use / G ff, the effective efficiency of 
energy injection, whereas for decays we define an alter- 
native quantity ?y. We then detail how likelihoods and 
corresponding constraints are obtained as a function of 
/eg or rj, including a simple prescription for combining 
results from multiple channels. 

Here we explicitly build upon existing public datasets 
[14, 17, 24], adding the necessary final steps for earlier 
results to be immediately implemented in, for example, 
multi-messenger dark matter analyses [26] and global fits 
to beyond the Standard Model particle theories [27-29]. 
This complements and extends recent efforts to construct 
public likelihood functions from Fermi-L AT [30], HESS 
[31] and IccCube [32] searches for dark matter annihila- 
tion. 

Annihilating dark matter. We will first consider 
the case of dark matter annihilating into standard model 
particles. Bounds on annihilation cross sections can be 
encoded in an integral involving f(z) and a set of prin- 
cipal component basis functions ei(z) that should be op- 
timized for annihilations, and which depend somewhat 
upon which experiment is being considered since they 
arc eigenvectors of the Fisher matrix. As discussed in 
[14, 17], f(z) can be expanded as 

oo 

ef{z)=Y,e i e i {z) (3) 
»=i 

where e = (av)/m x and Si = e f(z)- ei(z)/ei(z)- Ci(z). 
The inner product is the integral over z using the inte- 
gration limits Z\ = 86.83, z 2 = 1258.2. For annihilating 
dark matter, these basis functions arc chosen to maximize 
sensitivity to a general expected z-dependence for energy 
injection from annihilating dark matter, (based upon the 
energy density (1 + z) 3 S!dmPcC 2 of the dark matter itself) 
in the sense that the most important contributions to the 
observable energy deposition are described by the lowest 
components. The first three e^(z) are plotted in fig. 1. 

In fact, the contribution of the first component has 
been demonstrated to dominate, especially in the case of 
WMAP data, in the computation of the likelihood func- 
tion that is relevant for constraining annihilating dark 
matter models. Under the approximation of a Gaussian 
likelihood, and that the response of the CMB to the en- 
ergy injection is linear, the chi-squared is given by [14, 17] 

A x 2 = ]TA^/e- 2 (4) 

i 

where is the eigenvalue of the Fisher matrix corre- 
sponding to ei, for the relevant CMB experiment, and e 
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FIG. 2: 95% confidence level (CL) limits on cross section (av) for annihilation into the 15 channels indicated, for WMAP7 
(left, corresponding to (6)) and anticipated Planck results (right), for 14 months of observing time. As indicated by the 
numerals, the limiting curves fall into nearly-coinciding groups (1,4), 11, (7,8,9,12,10,15), (13,14), (2,5), (3,6), according to 
the channel-numbering key indicated on the figures. Horizontal line shows the nominal thermal relic density cross section, 
(av) = 3 x 10 _26 cm 3 /s. 



is a fiducial value, taken to be 2 x 10 27 cm 3 s 1 GeV 1 . 
For WMAP7, Ai = 0.279, while for Planck, A, = 
3.16, 0.691, 0.162 . . . , and the inclusion of contributions 
from the first three principal components is justified. 
Then, for example, the 2a upper limit on the cross sec- 
tion is given by 



(av) < 



2m x s 



VEiMCf-eO/fo-ei))* 



(5) 



The numerical value of the numerator in (5) (as well as 
that of (6) below) is correlated with the choice of nor- 
malization of the ei(x)'s. We use the functions shown in 
fig. 1, which are normalized such that [e.% ■ ei) = 23.9. 

However, it was observed in ref. [14] that the approx- 
imation of linear response is not yet very accurate for 
the magnitudes of £j that are allowed by the WMAP7 
data; instead one should perform a full likelihood analy- 
sis using CosmoMC [33] (and it is found that only the first 
principal component is required). This gives a somewhat 
weaker constraint than (5); the 95% c.l. upper bound is 



(av) < 1.2 x 10 



-26 



s-GeV (/ ■ ei)/(ei ■ ei) 



(6) 



(Throughout this paper, we use "95% c.l." to mean 2a, 
which is in fact 95.45% c.l..) 

In either (5) or (6), we see that the constraints can be 
expressed in terms of a single number involving integrals 
of f(z). To make contact with the earlier literature [10], 
we find it useful to consider a quantity f e s defined in 
terms of a "universal WIMP annihilation" curve e w (z) 
[14], / c ff ~ (/ • e w )/(e w • e w ), which has the interpreta- 
tion that / c ff < 1 denotes the average efficiency of energy 
injection for the annihilation channel of interest. It can 
be compared to the / mca n values tabulated in ref. [10] 



(column 2, table I). Here instead of defining / c ff directly 
in terms of e w , we use the expansion e w = J2i c i e ii an d 
the observation that using just the first term in the ex- 
pansion gives a good approximation. Thus for WMAP7 
we define 



fett 



(f-ei) 



(WMAP) 



(7) 



ci(ei • ei) 

where numerically c\ = 4.64. Then (6) can be expressed 



(av, 



< 0.26 x 10" 26 cm 3 s" 1 GeV _1 m x // ( 



off- 



(8) 



The analogous definition for Planck, which makes use of 
the contributions from the first three principal compo- 
nents, is 



1 



off 



/ATci 



1/2 



(Planck) (9) 



where we must use the Aj and ej appropriate to Planck. 
Then the 2a projected constraint from Planck takes the 
form (av) < 0.48 x 10 -27 cm 3 s -1 GcV _1 TO x // c ff . For an- 
nihilation into several channels i with branching fractions 
rj, the total / e ff is 



/off = ^2nf, 



eSA- 



(10) 



We present the 95% CL constraints as a function of 
m x for a range of annihilation channels in fig. 2, inter- 
polating from table I. The final states are pairs of e, [i, 
3, 3 V — > jj,, V — > t, q (light quarks: u, d or 



V 



3 V — > X denotes \X ~^ W followed by V — > XX with my 
1 GeV for X = e or fi and my = 5 GeV for X = r. 



4 



s), c, b, t, W, Z, g (gluons), 7, or h (Higgs). The lim- 
its from WMAP7 and projected constraints from Planck 
based on 14 months of observation time [34] are shown 
in the figure. Our 95% C.L. WMAP7 limits perfectly 
match those for e and /x found in present benchmark re- 
sults (e.g. at m x = 100 GcV, (<rv) < 4 x 10" 



for e + e 



{(TV 



< 1 x 10" 24 cur's" 1 for /Li+/i~ [13, 14]). 



We see approximately a factor of 6 improvement between 
WMAP7 and Planck, in agreement with the results of [12] 
but slightly less than the factor of 8 reported by [13]. 

Limits are strongest on the channels [e, V —> e] (which 
in fact coincide) and 7, which also coincides with the for- 
mer at masses > 100 GeV. The 77 final state limit has a 
different shape from the others because of the monochro- 
matic nature of the photons and the complex dependence 
of the efficiency of photon energy deposition on redshift 
and injection energy. This is encoded in the transfer 
function T 7 , and can be seen for example in fig. 1 of 
[17] and fig. 2 of [4]. At rcdshifts of z ~ 600, where 
the first principal component of fig. 1 is peaked, the 
Universe is approximately transparent to photons with 
energies E < 50 GeV, but rapidly transitions to com- 
pletely opaque for E > 100 GeV, due to scattering and 
pair production on CMB photons [4]. This causes the 
77 curve to track the electron curve above 100 GeV, as 
in both cases essentially all of the annihilation energy 
goes into heating the primordial gas. The 77 channel is 
the only one in which photons dominate over electrons 
in the final state spectra, and the only one to produce 
monochromatic photons - so it is the only one to show 
this complicated behaviour in the limits. 

The next strongest constraints are on the channels in- 
volving quarks, gluons and Higgs. These all coincide with 
each other (except that the top quark and Higgs have a 
higher mass threshold). The next weakest constraints are 
on [W, Z], followed by V — s- jj] and finally [r, V -¥ r\. 
These tendencies are indicative of the greater fraction of 
energy ending up in neutrinos for the least constrained 
channels. 



We stress that our results do not depend in any intrin- 
sic way upon the somewhat arbitrary definitions (7,9). 
Any quantity proportional to / • e% in the case of WMAP 
or (£,1 V / Aj(/ • e,/ej • e^) 2 ) 1 / 2 in the case of Planck would 
suffice to encode the necessary information for recover- 
ing the constraints. To this end, we list f e g values for 
both WMAP and Planck in table I, for the annihila- 
tion channels mentioned above and for the WIMP masses 
m x = 10, 30, 100, 300, 1000 GeV. They are in reasonable 
agreement with the / moan values tabulated in rcf. [10] in 
the cases that overlap. 



To determine constraints at an arbitrary confidence 
level, one would like to know the likelihood function for 
the annihilation cross section, assuming branching frac- 
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FIG. 3: First three principal components for decaying dark 
matter, from http://nebel.rc.fas.harvard.edu/epsilon/. Solid 
curves are for WMAP 7 and dashed for anticipated Planck 
data. 



tions Ti to channel i. For Planck, this is given by 
]nC((av)\m x ,ri) = 

-^/e 2 ff (m x ,n)A lCl ( 2xl0 ^ m 3 s -x ) (^) 



(11) 



from (5) and (9). As mentioned above, this expression, 
which assumes linear response of the CMB to the de- 
posited energy, is not accurate for WMAP. This can 
be corrected by making the replacement 2 x 10 -27 — > 
3.2 x 10 -27 in (11) (and one must also use the appro- 
priate value Ai = 0.279 for WMAP). Similarly, caution 
should be exercised if (11) is used to compute WMAP 
likelihoods far from 95% c.l. as the Gaussian approxi- 
mation was observed to be poor for WMAP likelihoods 
in ref. [14]. The WMAP 7- year limits can also be im- 
proved by ~15% if ACT data are added to the CosmoMC 
fit [13], and presumably slightly more than this if SPT 
data arc added. 

Decaying dark matter. A similar procedure can be 
used to constrain the fractional mass abundance 5Q of 
some metastable species present at the time of recom- 
bination, which decays into the same set of final states 
as we assumed above for annihilations. If the decaying 
species originally contributes Qj to the mass density of 
the universe relative to p c , and if the decay of a particle 
of that species converts a fraction <\> of its rest mass to 
standard model particles (and a fraction 1 — <f> to some 
lighter dark species), then <5f2 = cj&li. The injected elec- 
tromagnetic power density goes as 



dE 
dtdV 



l^-e- 1 ^ 5n Pc c 2 {l + zf. (12) 



The factor of (1 + z) 3 ~ n x as opposed to (1 + z) 6 ~ n 2 
for annihilations is the reason for the different powers 
of (1 + z) appearing in footnote 2. Unlike the case of 
annihilation, where the cross section appears only as a 
prefactor, here we have dependence upon the lifetime r 
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FIG. 4: 95% c.l. constraints on 5£1/£Iom versus r for x^> ee 
decays, with m x = 10, 100, 1000 GeV (solid, dashed, dotted 
lines respectively), from WMAP7 (upper curves) or Planck 
(lower curves). 

not only in the prefactor, but also in the z-dependent 
function e~*( z '/ T , which is not present in (1). This makes 
the analysis of decays more cumbersome than that of 
annihilations, because the transfer functions Ti(z' \z,E) 
now also depend upon r. We therefore need to compute 
(2) not only for all masses and final states, but for a grid 
of lifetimes as well. 

Note that this form of analysis allows us to not only 
constrain the lifetime and abundance of decaying dark 
matter candidates, but also the initial abundance of other 
metastable species, with such short lifetimes that they 
would not contribute to the current-day abundance of 
dark matter. For lifetimes less than ~10 17 -10 18 s the 
particles we refer to are therefore not dark matter in the 
usual sense, but rather some generic metastable particles. 

In the approximation of linear response to the injected 
energy, the chi-squared is most conveniently written in 
the form [17] 

Ax 2 =(/UV, (13) 

where we define 




to be the single number that encapsulates the energy 
injection history for a given mass, lifetime and final state 
(as / e ff does for annihilation). The likelihood is then 
simply 

ln£(r|m X! ri) = --(- ) t] 2 {t, m x , n). (15) 

The principal components e^(z) in rj are now special- 
ized to the case of decays, as well as depending upon 
the experiment. In contrast to annihilations, the first 



three components are needed for good accuracy, even for 
WMAP. We plot them for both experiments in fig. 3. The 
variances are given by at = (5.3, 6.4, 9.5 . . . ) x 10~ 25 s -1 
for WMAP7 and (1.2, 1.6, 2.3 . . . ) x 1CT 25 s" 1 for Planck 
[17]. The limits of integration appearing in / • a are 
z\ = 10, z<x = 1258.2, because n x redshifts slower than 
n 2 , , so lower redshifts are needed than in the case of an- 
nihilations. The inner products (and therefore the limits 
of integration) are now defined in terms of integrals over 
lnz rather than z, with normalization (e^ • e{) = 0.0543. 

Like /eff, V can be easily obtained for arbitrary branch- 
ing fractions r*j into multiple final states by taking 
the appropriate linear combination of values for single 
channels r\i 

v = yVjTft. (i6) 

i 

We have computed rj for each of the 15 final states 
that were considered as annihilation channels in the pre- 
vious section, over a grid of metastable particle masses 
10, 30, 100, 300, 1000 GeV as before (we augment this se- 
lection in the 7, t, W, Z and h channels where it would 
otherwise be too sparse), and a grid of 15 lifetimes. The 
results are given in table II. For the b channel, our results 
extend only as low as m x = 12 GeV, as all yields of [24] 
for the b channel go to zero as the energies of the the de- 
cay products approach 5 GeV (presumably because the b 
mass was approximated to 5 GeV in the calculations of 
[24]). 

Our lifetime grid is sufficiently dense to show the struc- 
ture of the limit on SD, as a function of r. For longer life- 
times, rj is the same as the values we give for r = 10 17 , 
allowing constraints and likelihoods to be extended to 
arbitrarily long lifetimes. An example of the resulting 
constraints is given in fig. 4 for the x ~ * e & channel. Our 
results for other channels agree with those presented by 
Slatyer [17] (as expected, given that we adopt her trans- 
fer functions here). In fig. 5 we plot the lower limit on 
the lifetime of a dark matter species decaying into any of 
the channels under consideration, assuming that it con- 
stitutes the entire presently-observed DM density. 

We again added extra points to the mass grid for the 7 
channel, in order to resolve the somewhat sharp feature 
seen in the region 100 < m x < 300 GeV, as discussed pre- 
viously in the results for annihilating dark matter. Note 
here that the mass range of the transition is approxi- 
mately a factor of two larger than seen in annihilation 
(as decays produce photons with energies equal to half 
the DM mass), but that this factor is not exact because 
the redshifts of importance differ in the two cases, so the 
energy above which the Universe is opaque to photons 
injected at relevant redshifts differs slightly. 

Our results for decaying species are based on the 
Fisher-matrix analysis of [14, 17], which assumes that 
the linear approximation for the response of the CMB to 
energy injection at high rcdshift holds. Given the small 
ionization fraction at the redshifts most important for de- 
cay (at least when the lifetime is long) we do not expect 
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FIG. 5: 95% c.l. lower bounds on the lifetime of decaying DM that accounts for the total DM density, from WMAP7 (left) 
or projected Planck data (right). The numbering of the final state channels is the same as in fig. 3: l:e, 2:/i, 3:r, 4:(V — > e), 
5:(V -> fj), 6:(V -> r), 7:(q = u, d, s), 8:c, 9:b, 10:t, 11:7, 12:5, 13:W, U:Z, 15:h. 



nonlinear corrections to be crucial. In principle though, 
this assumption needs to be verified, for example with 
explicit CosmoMC calculations, as was done for annihila- 
tion. 

It is also worth noting that the principal components 
ei(z) that we use for decay have been optimized for long 
decay lifetimes, leading to large leading contributions at 
relatively low redshift (fig. 3). In principle the ei(z) could 
also be optimized for different decay lifetimes, such that 
the optimal components for shorter lifetimes would prob- 
ably begin to resemble those for annihilation to some de- 
gree (fig. 1), potentially allowing the limits we give here 
to be improved for the shortest lifetimes. 

With tables I and II, and eqs. 11 and 15, we have pro- 



vided a fast and easy means for implementing CMB con- 
straints and likelihoods in future analyses of annihilating 
and decaying dark matter models, as well as some mod- 
els with metastablc particles that do not constitute dark 
matter. Our results can be interpolated to provide like- 
lihood functions and limits for arbitrary particle masses, 
and arbitrary annihilation and decay final state mixtures. 
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1.42 
1.43 
1.43 
1.44 
1.45 
1.69 
1.70 
1.72 
1.74 
1.75 
1.79 
1.80 
1.81 
1.81 
1.82 
1.84 
1.75 
1.76 
1.77 



0.86 
0.94 
1.02 
1.08 
1.09 
1.29 
1.33 
1.41 
1.47 
1.51 
1.42 
1.44 
1.50 
1.55 
1.56 
0.84 
0.90 
0.99 
1.05 
1.08 
1.29 
1.32 
1.39 
1.46 
1.52 
1.42 
1.43 
1.47 
1.53 
1.57 
1.25 
1.28 
1.31 
1.34 
1.38 
1.25 
1.27 
1.30 
1.34 
1.37 
1.24 
1.26 
1.29 
1.32 
1.36 
1.35 
1.36 
1.38 
1.39 
1.06 
1.11 
1.19 
1.04 
1.05 
1.07 
1.07 
1.08 
1.10 
1.25 
1.27 
1.30 
1.34 
1.37 
1.37 
1.39 
1.43 
1.39 
1.41 
1.46 
1.33 
1.36 
1.38 



0.53 
0.64 
0.75 
0.84 
0.85 
0.93 
1.01 
1.12 
1.21 
1.27 
1.10 
1.13 
1.21 
1.29 
1.32 
0.49 
0.59 
0.70 
0.80 
0.85 
0.92 
0.98 
1.08 
1.18 
1.27 
1.08 
1.11 
1.17 
1.25 
1.32 
0.94 
0.97 
1.01 
1.06 
1.11 
0.93 
0.96 
1.00 
1.05 
1.10 
0.91 
0.94 
0.98 
1.03 
1.09 
1.06 
1.07 
1.10 
1.11 
0.81 
0.89 
0.98 
0.78 
0.80 
0.82 
0.83 
0.84 
0.86 
0.93 
0.96 
1.00 
1.05 
1.10 
1.08 
1.11 
1.17 
1.10 
1.12 
1.19 
1.03 
1.07 
1.10 



0.29 
0.43 
0.56 
0.67 
0.69 
0.68 
0.78 
0.92 
1.03 
1.11 
0.86 
0.92 
1.01 
1.11 
1.15 
0.25 
0.36 
0.50 
0.62 
0.69 
0.67 
0.74 
0.87 
0.99 
1.10 
0.85 
0.88 
0.96 
1.06 
1.15 



0/ 

0.75 

0.79 

0.85 

0.92 

0.71 

0.74 

0.79 

0.85 

0.91 

0.68 

0.72 

0.76 

0.82 

0.90 

0.86 

0.87 

0.90 

0.92 

0.65 

0.74 

0.84 

0.60 

0.62 

0.64 

0.66 

0.67 

0.70 

0.71 

0.74 

0.79 

0.85 

0.91 

0.87 

0.91 

0.98 

0.89 

0.92 

1.00 

0.82 

0.87 

0.91 



0.13 
0.28 
0.44 
0.56 
0.59 
0.52 
0.62 
0.78 
0.91 
1.01 
0.71 
0.77 
0.88 
0.99 
1.04 
0.09 
0.21 
0.37 
0.50 
0.58 
0.50 
0.58 
0.72 
0.86 
0.98 
0.68 
0.73 
0.82 
0.93 
1.03 
0.55 
0.59 
0.65 
0.72 
0.79 



0.55 
0.59 
0.65 
0.71 
0.79 
0.51 
0.56 
0.62 
0.68 
0.77 
0.72 
0.73 
0.77 
0.78 
0.54 
0.64 
0.75 
0.48 
0.51 
0.53 
0.55 
0.56 
0.60 
0.55 
0.59 
0.65 
0.72 
0.79 
0.73 
0.77 
0.86 
0.74 
0.78 
0.87 
0.68 
0.74 
0.78 



0.02 
0.19 
0.36 
0.49 
0.53 
0.40 
0.52 
0.69 
0.83 
0.94 
0.60 
0.67 
0.79 
0.91 
0.97 
-.03 
0.11 
0.28 
0.43 
0.52 
0.38 
0.47 
0.62 
0.77 
0.91 
0.57 
0.63 
0.72 
0.85 
0.96 
0.44 
0.49 
0.55 
0.63 
0.70 
0.44 
0.49 
0.55 
0.62 
0.70 
0.40 
0.45 
0.52 
0.59 
0.68 
0.62 
0.64 
0.68 
0.70 
0.48 
0.58 
0.70 
0.40 
0.43 
0.46 
0.48 
0.49 
0.53 
0.44 
0.49 
0.55 
0.63 
0.70 
0.64 
0.68 
0.78 
0.65 
0.69 
0.79 
0.58 
0.65 
0.69 



-.05 
0.12 
0.30 
0.45 
0.49 
0.33 
0.45 
0.63 
0.78 
0.89 
0.53 
0.60 
0.73 
0.86 
0.93 
-.10 
0.04 
0.22 
0.38 
0.48 
0.30 
0.39 
0.55 
0.71 
0.86 
0.49 
0.56 
0.66 
0.79 
0.91 
0.37 
0.42 
0.49 
0.56 
0.65 
0.36 
0.42 
0.48 
0.56 
0.65 
0.32 
0.38 
0.45 
0.52 
0.62 
0.56 
0.58 
0.62 
0.64 
0.44 
0.54 
0.66 
0.35 
0.38 
0.41 
0.43 
0.45 
0.49 
0.37 
0.42 
0.49 
0.57 
0.65 
0.57 
0.62 
0.72 
0.58 
0.62 
0.73 
0.52 
0.59 
0.63 



-.09 
0.08 
0.27 
0.42 
0.47 
0.28 
0.40 
0.59 
0.75 
0.87 
0.48 
0.56 
0.69 
0.83 
0.90 
-.15 
-.00 
0.18 
0.34 
0.45 
0.25 
0.35 
0.51 
0.68 
0.84 
0.44 
0.51 
0.62 
0.75 
0.88 
0.31 
0.37 
0.44 
0.52 
0.61 
0.31 
0.37 
0.44 
0.52 
0.61 
0.26 
0.33 
0.40 
0.48 
0.59 
0.52 
0.54 
0.58 
0.60 
0.41 
0.52 
0.64 
0.32 
0.35 
0.38 
0.40 
0.42 
0.47 
0.32 
0.37 
0.44 
0.53 
0.61 
0.53 
0.58 
0.69 
0.54 
0.58 
0.70 
0.47 
0.55 
0.60 



TABLE II: Values of - log 10 r? = - log 10 E i ((/ • eO/M^ • e t )f} 1/2 (where is in units of 10 2o s 1 ) as a function of log 10 r 
(with r in seconds) and m x in GeV, for each of the 15 decay channels whose numbering corresponds to the key in figs. 2 or 5. 
Figures on left for WMAP7 and on right for Planck. Values for longer lifetimes are the same as for log 10 r = 17. 
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